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Abstract. Within the formalism of the matrix product ansatz, we study a two- 
species asymmetric exclusion process with backward and forward site-ordered se- 
quential updates. This model describes a two-way traffic flow with a dynamic im- 
purity and shows a phase transition between the free flow and the traffic jam. We 
investigate characteristics of this jamming and examine similarities and differences 
between our results and those with the random sequential update [lj . 

1 Introduction 

A variety of phenomena can be modeled by the one dimensional asymmet- 
ric simple exclusion process (ASEP) and its generalizations (see &|| and 
references therein). The model has a natural interpretation as a description 
of traffic flow and constitutes a basis for more realistic ones |§|. In traffic 
flow theories, the formation of traffic jams due to "impurities" is one of the 
fundamental problems. 

In the ASEP models, two kinds of impurities are discussed in the liter- 
ature. The first one is "dynamic" impurities, i.e., defective particles which 
jump with a rate lower than others In the traffic terminology, such 

moving defects can be visualized as slow cars on a road. The other kind is 
"static" impurities such as imperfect links where the hopping rate is lower 
than in other links [p pp| , pl] , p2| . Both types of impurities can produce shocks 
in a system. Presently a limited amount of exact results is available for the 
shock formation and most of them are for models with the random sequential 
update |,|@. 

The implementation of the update is an essential part of the definition of 
a model and it is of prime interest to determine whether distinct updating 
schemes can produce different behaviors. The aim of this work is to investigate 
consequences of changing the updating scheme of the model. Here, we study 
an exactly solvable traffic model with two types of ordered sequential updates, 
which is identical to the model studied in Ref. [Q except for the updating 
schemes. 
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2 Model Definitions and Matrix Product Ansatz 
2.1 Two- Way Traffic Model 

Consider two parallel one dimensional chains, each with N sites and the 
periodic boundary condition. There are M cars and K trucks in the first and 
the second chain, respectively, and cars move to the right while trucks move 
to the left. We introduce inter-chain interaction that forbids a car and a truck 
to occupy two parallel sites simultaneously. Then the state of the system can 
be described by a single set of occupation numbers (ti, T2, • • • , tjv) where 
Ti = (empty site), 1 (occupied by a car), or 2 (occupied by a truck). There 
are three possible hopping processes: 



The reduction factor (3 (1 < (3 < oo) is related to the width of roads: (3 = 1 
corresponds to a very wide road or a highway with a lane divider, and (3 = oo 
corresponds to a one lane road. Thus the value 1 — 1/(3 can be used as a 
measure of the road narrowness. 

Recently the model is studied with the random sequential update (RSU). 
In an infinitesimal time interval dt, each link is updated with the probability 
that is the product of the relevant rate and dt [Since rj does not affect the 
dynamics at all in the RSU scheme, we choose rj = 1 for the RSU scheme]. 

Alternatively the ordered sequential update (OSU) can be used. One 
first chooses a particular site, i.e. the site N, and updates the state of the 
links consecutively either in the backward direction (N, N — 1), (N — 1,N — 
2), ■ • • , (1, N) [backward sequential update (BSU)] or in the forward direction 
(N, 1), (1, 2), • • • , (N - 1, N) [forward sequential update (FSU)]. In contrast 
to the RSU update, the time is discrete in the BSU and FSU schemes, and in 
each time step, hopping occurs in links with the probabilities that are equal 
to the rates in Eq. (|l|). 

2.2 Matrix product state 

The two-way traffic model is equivalent to a two-species ASEP and it can be 
solved exactly by the method of the matrix product state (MPS). In the RSU 
scheme, the steady state weight P s of a given configuration (ti,T2, • • • , tjv) is 
proportional to the trace of the normal product of some matrices [0 : 



where Xi = D for Tj = 1; E for rj = 2; and A for Tj = 0, and these matrices 
satisfy the quadratic algebra 



(i) Car hopping with rate rj : (1, 0) — »■ (0, 1) 

(ii) Truck hopping with rate 777 : (0, 2) — ► (2, 0) 

(iii) Car-truck exchange with rate 5 : (1, 2) — > (2, 1) 



(1) 



P s (n,T 2 , ■ ■ ■ ,t n ) ~ Tr(X!X 2 ■■■X N ) 



(2) 



DE = D + E, aAE = A, (3D A =A {a = (3-y) . 



(3) 
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In a similar way, P s in the BSU (FSU) scheme can be written as follows: 

P s {n,T 2 , ■ ■ ■ ,rv) ~ Tr(XiX 2 • • -X N ) (4) 

where Xi — D, E, or A (Xn = D, E, or A) depending on T} (tjv). Note that 
the matrices at the site N are different from those at other sites jL3|. The 
presence of the hatted matrix in Eq. (|j) breaks the translational invariance 
of the problem. In Ref. ]l4|] , it is shown that with a simple assumption, 

A = A + a, D = D + d, E = E + e, (5) 

where a, d, and e are carefully chosen real numbers, the relevant matrix al- 
gebra for A, D, E becomes identical to the one in Eq. (|^) except for the 
renormalizations of a and to a and [3. In the BSU scheme, the choice 
a = 0, d = — T]/0, e = r;/ (/3 — fy) gives a — a((3 — r?)/(/3 — a??), (3-/3, and in 
the FSU scheme, the choice a = 0, c? = r)/(/3 — rf), e = —r)/j3 results in a = a, 
p= {fl- r?)/(l -77). 



3 Average Velocities 

In this section we consider the special case where there is only one truck 
and evaluate average velocities. The MPS method allows exact evaluations 
of velocities. We first present the car velocity in the thermodynamic limit. In 
the BSU scheme, 

V -(l-n)if n/3<l 



1 — r/n 

(«car>= < , (6) 



and in the FSU scheme, 



77 1 — n 
— r) n 

1 — 77(1 — n) 



if n/3 > 1 



if n(3 < 1 



^(i^) if n ^>i 



(7) 



Figs. 1(a) and (b) show the behaviors of (v cal ) in the two updating schemes 
as a function of the road narrowness 1 — 4 while the values of 77, j(= ^) and 
n are fixed. In both schemes, a continuous phase transition is evident: While 
(fear) is constant below the critical narrowness, it begins to drop suddenly at 
the critical narrowness (n/3 = 1), generating a cusp. Thus above the critical 
narrowness, a single truck results in global effects. The functional dependence 
of the decrease differs in the two updating schemes. 

The exact expressions for the truck velocity, on the other hand, are rather 
lengthy especially below the critical narrowness and we present the functional 
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(d) 



Fig. 1. (a) Average velocity of cars in the BSU and RSU schemes for n = 0.3 
and different values of r\. (b) Average velocity of cars in the FSU and RSU 
schemes for n — 0.3 and different values of r\. (c) Average velocity of the 
truck in the BSU and RSU schemes for n = 0.3 and different values of rj. (d) 
Average velocity of the truck in the FSU and RSU schemes for n = 0.3 and 
different values of 77. After Ref. H]. 
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dependence on the road narrowness only through figures [Figs. 1(c) and (d)]. 
For exact expressions, see Ref. The appearance of a continuous transition 
at the critical narrowness {n(3 = 1) is clear and the functional dependence of 
the truck velocity again differs in the two updating schemes. 

4 Density Profile 

In the RSU scheme, the probability to find a car at a site depends only on its 
relative distance x to the truck due to the translational invariance. In Ref. §, 
it is found that in the high density phase n/3 > 1, the probability or density 
profile (n(x)) becomes 



Note that the system consists of two regions, a traffic jam region in front of 
the truck and a free flow region behind it. In the low density phase n/3 < 1, 
on the other hand, the presence of the truck has only local effects. In the 
thermodynamic limit, the car density becomes 



which shows that the disturbance by the truck decays exponentially with a 
characteristic length scale | ln(n/3)| _1 . 

In the ordered sequential updates, complications occur in the definition of 
the probability itself. Since the choice of a particular site as a starting point 
of the update breaks the translational invariance, the probability to find a car 
at x sites in front of the truck depends not only on the relative distance x but 
also on the truck location. This unnecessary complication can be avoided by 
choosing the starting point in an even way. In Ref. []l4| , it is found that with 
this "homogeneous" choice, the expressions for the probability in the BSU 
and the FSU scheme become identical to Eqs. except for the replacement 
of the bare parameters a and (3 with the renormalized ones a and j3. 

5 Concluding remarks 

Characteristics of an exactly solvable two-way traffic model are investigated 
with the OSU scheme and both qualitative and quantitative differences are 
found compared to those with the RSU scheme |lj . Our approach is based on 
the so-called matrix product formalism which allows analytic solutions. In the 
OSU schemes, the choice of a particular site as a starting point of the update 
breaks the translational invariance of the steady state measure, which is also 





(9) 
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evident in the form of the MPS. Thus an averaging over the different choices of 
the update starting point is necessary to restore the translational invariance 
to the system. Performing the translationally invariant averaging, it is found 
that for characteristics such as the density profile of cars and the density- 
density correlation function, the differences in the updating schemes can be 
fully taken into account by the proper renormalization of the parameters 
a and (3. However this is not the case with average velocities. Changing 
the update scheme affects velocities in a more complicated manner and the 
renormalization of the parameters is not sufficient to account for different 
behaviors of (i> car ) and (utruck) in different updating schemes. Especially the 
dependence of (v cax ) and («truck) on the road narrowness 1 — can vary 
qualitatively depending on the updating schemes. 
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